arXiv:1508.01399v2 [cond-mat.str-el] 10 Dec 2015 


Local realizations of anyon exchange symmetries without lattice 

dislocations 

Miguel Jorge Bernabe Ferreira,®* Pramod Padmanabhan,®’^! 

Paulo Teotonio-Sobrinho®! 

December 14, 2015 


“ Departamento de Fi'sica Matematica, Universidade de Sao Paulo 
Rua do Matao Travessa R Nr.187, Sao Paulo, CEP 05508-090 
^ Department of Physics and Astronomy, University of California, Irvine 
4129 Frederick Reines Hall Irvine, CA 92697-4575 

Abstract 

The global e-m exchange symmetry of the toric code is realized locally through an exactly solvable 
Hamiltonian on a two dimensional lattice which has no lattice dislocations and their associated defect line. 
The Hamiltonian is still changed locally in selected sites where we wish to realize this anyon symmetry. We 
refer to these selected sites as defect sites in analogy with the usual lattice defects. The operators on the 
defect sites condense dyons of the toric code and are shown to support states which obey the fusion rules of 
Ising anyons just as the lattice dislocations thus achieving the transition to the Ising phase from the toric 
code phase. They can also be introduced in an entire region leading to an idea of non-localized defects. 
The method leads to a natural generalization for other Abelian groups where they help realize all the anyon 
exchange symmetries locally. 


Introduction 


Physics in three dimensional space is governed by bosons and fermions as the fundamental group is the 
permutation group, 5'„ [T] which have two irreducible one dimensional representations. The story changes in 
two dimensions where the fundamental group is given by the braid group, of n particles which leads to a 
rich variety of statistics. These are termed as anyons [2] which have recently found applications in the held of 
topological quantum computation HIS]. Exactly solvable lattice models have been used to obtain these particles 
as low energy excitations of which the most famous one is the toric code model of Kitaev and it’s variants mn\- 
These are examples of systems with long-ranged entangled ground states in which information in qubits could 
be encoded non-locally. Soon they were realized to be a major part of the classihcation of phases of matter at 
low temperatures and were dehned to possess a new type of order called topological order [3]. 

The phase of the toric code model is described by it’s anyon content which are the representations of the 
quantum double of Z 2 [5]. They are the electric charge, e, magnetic hux, m, the combination of the two 
called the dyon, e and the vacuum, 1 which can be thought of as a trivial anyon. They obey the fusion rules 
exm = e; exe = mxm = €Xe = l; exe = m ; mxe = e with the other fusion rules being the trivial ones. The 
dyon is a real fermion and the rest are bosons. It is easy to see that these rules are invariant under an interchange 
of e and m labels. This is the electric-magnetic duality of the toric code model which was known earlier [S]. 
In recent literature this global exchange symmetry was seen to be part of the anyon symmetry group mm- 
Gauging this symmetry leads to new topologically ordered phases from a given parent phase with interesting 
features such as the emergence of non-Abelian anyons from Abelian phases. Recently efforts have been made to 
construct models realizing these symmetries locally with the earliest examples being the work of Bombin m 
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where a microscopic model with lattice defects were introduced to gauge the electric-magnetic symmetry of the 
Z 2 toric code. This process led to the emergence of Ising anyons in a model where the underlying degrees of 
freedom belonged to that of the toric code. Thus this served as an example of realizing non-Abelian anyons 
from an Abelian system |12j . Although Abelian systems are simpler to be experimentally implemented, they 
are not enough for quantum computation m Finding Abelian systems supporting non-Abelian anyons not 
only deepens our understanding of topological phases but is also relevant for quantum computation [TillTSHTB] . 
This model was further generalized for the Z„ case in m- Ising anyons were shown to be obtainable from the 
toric code anyons as a fusion category by superposing toric code anyons in [HI HI]. These features were also 
studied for the toric code in [19]. These efforts were taken further in a more complete theory of such defects 
in Ablelian theories in the effective field theory language in [50] where such symmetries were also realized in 
bi-layer fractional quantum Hall states |21] . A mathematical theory was developed recently under the name of 
G-crossed braided fusion categories in |5]. Recently global realizations of these symmetries were obtained as 
models for symmetry enriched topological (SET) phases in [521 [53]. Despite these important advances there 
are still very few exactly solvable microscopic lattice models that realize these phases locally in a systematic 
manner. 

In this article we propose a microscopic, exactly solvable two dimensional lattice model which just uses the 
familiar operators appearing in the description of the Z 2 toric code to realize the electric-magnetic symmetry 
locally. Moreover we also show that Ising anyons emerge in this phase. We do not introduce any physical 
defect/dislocation in the lattice to achieve this. However we change the Hamiltonian at sites, s = (v,p) where 
the electric-magnetic symmetry is realized locally. We call these sites as defect sites and they should not be 
confused with the lattice defects used in m- The operators associated with these defect sites were introduced 
earlier in [23]. These defect sites act as sources and sinks for the dyons. The eigenvalues of the operators 
associated to the defect sites define two sectors just as in the case of the defects of m- The difference in this 
model is that the defect site achieves both the exchange of e and m and when we have two such defect sites we 
can realize the Ising anyon phase as well. This is in contrast to the lattice defects introduced in m where the 
defect operators at the ends of the branch cut are very important for realizing the Ising anyon phase. Another 
important difference is that we can now create defect regions which consists of an area of defect sites. This helps 
switch the charge and flux in a non-local way. When we have two such defect regions then we can fuse them in 
the same way as the localized defect sites to again realize the Ising anyon phase. By comparing the degeneracies 
of this model with those using lattice dislocations [13] we obtain a relation between the two approaches. Finally 
the model presented here can be easily generalized to other Abelian cases and possibly non-Abelian ones as 
well. In particular the model naturally realizes all the anyon exchange symmetries of the Zjv toric code anyons. 


The model 

The setup of the model is similar to that of the toric code with the degrees of freedom located on the links of 
a two dimensional lattice. For simplicity we choose this to be a square lattice though the construction goes 
through for a lattice with an arbitrary triangulation just like the toric code model. 

The dynamics is given by vertex Ay and plaquette Bp operators, defined as the usual toric code stabilizers 

Ay = af and Bp = (^ af . (1) 

JGStar(i;) l&9p 

Apart from these operators we also consider operators that act on the sites of the lattice, called site operators 
Os- Each site s = (u,p) is formed by a vertex and a plaquette as shown in figure [TJ The site operator is define 
as 

Os' = Ay)Bp) , where s' = {v',p') . 

The action of this operator on the qubits is shown in figure [T] 

Consider S = {sq, si, • • • , s„}, a set of sites of the lattice. The Hamiltonian is given by 

= ( 2 ) 

sGS P^S v^S 

where hy p ^ S {v ^ S) we mean the plaquette p (vertex v) is not part of any site which belongs to S. It is easy 
to check that this model is exactly solvable and it’s spectrum can be readily written down just as in the toric 


2 



Figure 1: Action of the operator Og on a site s = (v,p). The cr® operator acts on the spikes of the vertex v. 
The remaining Pauli matrices act on the sides of the plaquette p. 


code case. The ground states satisfy Note that the ground states of the 

toric code are also ground states of this model. The difference arises due to = |^'°), which can also be 

written as A„i3p|d'°) = |4'°) (s = {v,p)) which implies or = i?p|^'°) = —|1E''’). 

This implies that the dyon on site s € S gets condensed into the ground state. On a closed manifold the 
toric code excitations can only be created in pairs, thus if a given state has an even number of dyonic 
excitations localized on the sites s G S (and no other kind of excitation anywhere else) it is still in the vacuum 
of the Hamiltonian ©• The ground state degeneracy of this model is given by the sum of all possible ways 
of creating a pair of dyons on the n sites of S and is given by 2^3+”“^. In the case of the torus this reduces 
to 2"+^ for a system with n defect site operators. In earlier realizations of the electric-magnetic duality using 
lattice dislocations the ground state degeneracy was found to be 2“+^ for m lattice dislocations [13] and this 
also implies that the dislocations carry non-integral quantum dimension hinting at their non-Abelian behavior. 
Comparing the two ground state degeneracies we see that n defect site operators correspond to 2n lattice 
dislocations which implies that each defect site operator mimics the behavior of two lattice dislocations. Since 
the defect sites are a combination of two lattice dislocations, we may ask how the phase with non-Abelian 
behavior arises in this model. We shall see that the model supports states of the form e -I- m apart from the 
usual toric code anyons which behave like Ising anyons as shown in m- 


Excitations 


The excitations of this model are similar to the excitations of the toric code except in the region where we 
have the defect sites that is the region where the Hamiltonian is given by the operators Og^ acting on the sites 
Si- In the rest of the lattice the excitations are the deconfined charges, e, fluxes, m and the dyons, e. In the 
Hamiltonian in ([2|) we have n sites Si where the operator Og^ acts. We will consider different cases starting with 
the case where we have a single isolated defect site, that is n = 1. Then we will consider the case with n > 1 
acting on isolated sites. Finally we will consider the case where we still have n > I, but now the defect sites are 
not isolated but grouped into a region covering an area A. In what follows we will only see the behavior of the 
excitations near the defect site or the defect region. 

Single defect site (n = 1) 


The site operator Og has two eigenvalues, ±1 with the corresponding states denoted by a±. They split the 
Hilbert space into two sectors with the excitations of the system behaving differently depending on which sector 
the state is in. 

Consider the creation of a sector where the eigenvalue of the site operator is —1 or the (t_ state. A way 
of doing this is by creating an isolated charge near the defect site which also excites the Og operator into the 
(T_ sector as shown in hgure 2(a) A similar situation occurs for the flux excitations of the system as shown in 
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(a) 



(b) 


Figure 2: An isolated charge or flux for which the a- sector is a source or sink of. The red strings correspond 
to string operators formed out of cr^ operators and the blue strings correspond to string operators made of 
operators. 


figure 2(b) The U- sector acts as a sink or source for the charge (and also the flux) as shown in figures 2(a) 


and 2(b) respectively. Thus in this sector the system can support a state of the form e - 
charge and a flux m- 


m that is an isolated 


The tT+ sector acts as source/sink for the dyonic excitations as shown in figure 3(a) Apart from acting as a 



Figure 3: The anyons that can occur in the (t+ sector, (a) an isolated dyon and (b) a pair of charges, and (c) 
the permutation of the e-m excitations. 


source/sink for the dyonic excitation, the (Tj^ sector can also host a pair of charges as shown in hgure [T(b)| or a 
pair of fluxes which are part of the deconfined excitations of the toric code. Thus in this sector the system can 
host a state of the form 1 + e. 

Another interesting feature of the a+ sector is that it can interchange a charge and a flux, as shown in 
figure 3(c) This can be seen in the following way. The charge and the flux excitation that compose the dyonic 
excitation shown in figure [3^a)] can be moved apart from each other giving rise to the state shown in figure 3(c 


This implies that in this sector of the Hilbert space a charge excitation can be made into a flux excitation by 
crossing the defect site. Note that this occurs without the introduction of any lattice dislocation. 


Isolated defect sites (n > 1) 


As in the case with n = 1 the defect sites act as sources and sinks for the toric code excitations according 
to the eigenvalues of the operators defined on those defect sites and as before each of them also switch a charge 
into a flux and vice versa. However the phase in which this model exists is different as the ground degeneracy is 
different from that of the toric code case unlike the case with a single defect site which has the same degeneracy 
as the toric code. The only new feature now is that we can “fuse” two defect sites by fusing the particles for 
which these defect sites are sources. 

The way to fuse the two sectors <j± is by moving the anyons for which these sectors are a source of into a 
common site of the lattice. We can call these anyons as defect excitations. To move a defect excitation one has 
to extend the string that ends at the defect site. Figure 4(a) shows a possible way of fusing (t+ and tT_ where 
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the respective anyons are moved into the shaded region. In the same way the U- excitation could have also 





Figure 4: Fusion rules of the (t+ and (t_ sectors, with the fusion channels of (t+ and (t_ giving (a) charge and 
(b) flux, and the fusion channels of a- and cr_ giving (c) a dyon and (d) the vacuum. 


moved into the shaded region by applying a blue string as shown in figure |4(b)[ in this case fusion would give 
e X TO = e, a charge excitation as the result. In a similar way we could have fused the vacuum in the (t+ sectors 
with the e or TO particles from the (t_ sector resulting in a e and to particle respectively. Thus there are several 
channels for the fusion of the CT-i- and (t_ sectors leading to the non-Abelian fusion rule 


a± X = e + m . 


Using the same procedure it is not difficult to see that the fusion rules for a± x are also non-Abelian and 
are given by 

cr± X cr± = 1 -I- e . 


Figures 4(c)| and |4(d)| illustrate these two ways for the fusion of two (7_. This comes from the fact that there 
are two ways to move U- (and also cr+) excitations since they are source/sink for two kind of particles as shown 
in figure [2] and |3(a)[ 

We can also fuse the anyons created far away from the defect site with the defects tT+ and tT_. By using the 
fusion rules of the toric code we easily obtain 


(T± X e = CT± ; cr± X e = (T± X to = (Tip. 

In particular we can choose cr = e -I- to and ip = e. Along with the vacuum I they give rise to the fusion 
rules of Ising anyons [Hi¬ 
lt is easy to generalize this statement for a system with n defect sites. The two sectors are given by the 
product of the eigenvalues of the operators Osi, • • • , Os„ on the defect sites, si, • • • , s„. 


Defect sites in a region A 


Finally we consider the situation where the defect sites are grouped into a region A with n sites. This model 
has the same ground state degeneracy as the case where the defect sites are isolated. The features of the defect 
sites are unchanged from the previous situations including the realization of the fusion rules of the Ising anyons. 
The two sectors in this case are given by the possible product of the eigenvalues of the operators Os ^, * G A 
which are ±1 as before. We can think of this whole region as a defect which is spread in a region with area A. 
Thus we have a non-localized defect unlike the previous two situations where we dealt with isolated, localized 
defect sites. Inside this region there is no distinction between e and to excitations. 

The new feature in this case which was not present in the earlier cases is that the transmutation of the 
charge and the flux can happen non-locally. By this we mean that a charge which enters a defect region A can 
exit it from any point of the region as a flux, to as shown in figure [5l That is there exists a string operator 
which helps us achieve this. Figured also shows an example of such a string operator. 

As before the two sectors are labeled by a+ and (t_ . If we have two isolated defect regions we can think of 
fusing them as we did in the case of isolated defects. We again have Ising anyons in this case as the argument 
for seeing this to be the case is similar to the previous situation. 
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Figure 5: Non-local transmutation of a charge, e and flux, m in a defect region A. The area A is represented 
by the shaded region. 


Finally we can think of the extreme case where we only have defect site operators all over the lattice. Then 
we are completely out of the Ising phase and are in a phase where there is no distinction between the charges 
and fluxes and all the dyons are condensed. 


Other ZjY groups 

An advantage of this approach is that they can be naturally generalized to the case of other groups and inputs, 
which the parent quantum double models are based on. In particular they can be easily written down for 
the Zjv case where we have N — 1 I .2 anyon exchange symmetries. The GSD of these models is given by 
]\[2g+J2i=i ni-i denotes the number of defect sites which realize the Ah anyon exchange symmetry. 

On the torus this reduces to the result in na, where n is a pair of dislocations realizing one of the anyon 

exchange symmetries. Each dislocation displays non-Abelian behavior while the defect site operators support 
superpositions of Z]\[ toric code anyons which have non-Abelian behavior. 


Remarks 

Decorating topological phases with dislocations were studied from the categorical point of view in m and 
their appearance can be understood as follows. Physical systems with boundaries have two phases, a bulk 
and a boundary phase. Depending on the type of the boundary, smooth or rough, we have different boundary 
phases. Domain walls between different topologically ordered phases contain a phase different from the phases 
on either side of it as the bulk anyons from the two phases fuse into the domain wall making them domain 
wall excitations. This problem reduces to the problem of identifying the boundary phase once we fold the 
system along the domain wall |19j . The notion of fusion of two such domain walls in the presence of three 
different topologically ordered phases into a single wall can then be introduced. Along similar lines joining two 
domain walls, between different topologically ordered phases, at a point results in a new phase locally at these 
points. These points have been studied as defect points or dislocations in the lattice. As mentioned before 
these are physical artifacts of the lattice on which the system lives on resulting in a change in the Hamiltonian 
locally around the regions where they are present. Thus the introduction of these lattice dislocations changed 
the Hamiltonian locally which described a new phase, in this case the Ising phase emerged from the toric 
code phase. The operators introduced in this paper achieve the effects of these defects/dislocations without a 
modification of the lattice through a new microscopic Hamiltonian which again described the Ising phase by 
condensing the dyons of the toric code phase. 
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